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ABSTRACT 


A  digital  computer  program  was  written  which  will  perform 
interior  ballistic  calculations  on  an  IBM  709  or  IBM  7090.  A  brief 
description  is  given  as  well  as  an  outline  of  the  method  of  analysis 
which  this  program  uses.  Examples  are  presented  and  both  input 
and  output  formats  are  discussed. 


CONCLUSIONS 

It  is  not  the  purpose  of  this  report  to  evaluate  or  verify  the  accuracy 
of  the  Hirschfelder  System  on  Interior  Ballistics.  Its  only  purpose  is 
to  describe  a  digital  computer  program  which  will  perform  the 
Hirschfelder  calculations. 

The  reader  will  note  that  the  results  of  the  hand  calculations  given 
in  Part  IV  correspond  exactly  to  the  computer  results  given  in 
Appendix  B. 


INTRODUCTION 


This  report  describes  a  digital  computer  program  for  performing 
the  interior  ballistic  calculations  of  J .  O.  Hir  schfelder . 

The  Hirschfelder  System  was  developed  between  1942  and  1945  and 
makes  these  basic  assumptions: 

1.  A  "first-degree  burning  law"  is  used  in  finding  solutions  to 
the  fundamental  ballistic  equations. 

2.  The  powder  gas  is  distributed  according  to  Kent's  solution 
of  the  problem  of  the  motion  of  the  powder  gas. 

3.  The  heat  lost  to  the  bore  up  to  any  instant  is  proportional 
to  the  square  of  the  velocity. 

4.  The  friction  of  the  projectile  is  taken  as  equivalent  to  a 
resisting  pressure  on  the  base  of  the  projectile  which  is  equal  to  a 
constant  fraction  of  the  average  pressure. 

This  report's  findings  are  divided  into  four  parts: 

The  first  part  discusses  the  input  and  output  formats. 

The  second  part  is  devoted  to  a  study  of  the  basic  equations 

used  in  the  program.  (It  is  assumed  that  the  calculations 

found  in  Part  II  will  be  used  in  conjunction  with  Reference  1.) 


The  third  part  is  a  brief  presentation  of  program  logic. 

In  Part  IV  an  example  is  given  in  which  Cases  2,  3  and  4 
are  presented. 

The  program  will  solve  four  different  types  of  problems  (cases). 
In  Case  1,  various  1  s  are  given  along  with  maximum  pressure  and 
the  computer  will  find  velocity  and  web.  In  Case  2,  charge  and  maxi¬ 
mum  pressure  are  given  and  web  and  velocity  are  found.  In  Case  3, 
charge  and  velocity  are  given  and  web  and  maximum  pressure  are 
found.  In  Case  4,  charge,  velocity  and  web  are  given  and  the  burning 
constant  is  found. 
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The  appendices  include  tables  giving  propellant  codes  and 
constants,  the  program  output  for  the  example  given  in  Part  IV  and 
also  the  complete  Fortran  Program. 

The  reference  used  is  a  revision  and  consolidation  of  seven  progress 
reports  on  interior  ballistics  written  by  J.  O.  Hirschfelder  and  others 
of  the  staff  of  the  Geophysical  Laboratory,  Carnegie  Institution  of 
Washington  (Reference  1), 
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DISCUSSION 


PART  I  -  INPUT  AND  OUTPUT  FORMATS 

Input  Format 


The  first  104  data'cards  will  be  the  same  for  any  group  of  runs. 
These  cards  incorporate  constant  tabular  values  (Part  III  of  this 
report) . 

Following  this  set  of  cards,  any  r  amber  of  runs  may  be  included 
provided  each  run. is  in  correct  sequence.  For  each  run,  the  follow¬ 
ing  format  is  used.  .  - 

Card  1 


Space  1  is  reserved  for  a  numerical  code  giving  the  type  of 
problem  to  be  solved. "1".  means  that  maximum  pressure  and  para¬ 
meters  4>'s  are  given  and  web  and  velocity  are  to  be  found.  "2" 
means  that  charge  and  pressure  are  given  and  web  and  velocity  are 
to  be  found.  "3"  means  that  charge  and  velocity  are  given  and 
maximum  pressure  and  web  are  to  be  found.  "4"  means  that  charge, 
velocity  and  web  are  given  and  the  burning  constant  is  to  be  found. 

Spaces  2  and. 3  are  reserved  for  a  numerical  code  which  denotes 
the  propellant  to  be  used  in  accordance  with  Table  1.  The  propellant 
constants  are  (given  in  Table  2)  automatically  selected  from  the 
memory  of  the  machine.  If  no  propellant  code  is  specified  or  if  it  is 
given  as  "99,"  these  constants  must  be  given  as  part  of  the  data, 
appearing  on  Card  2  below. 

Space  4  is  reserved  for  the  code  form  function.  A  "1"  denotes  a 
propellant  of  single  perforation.  A  "2"  denotes  a  multiperforated 
propellant  (for  example,  seven-perforated  grains). 

Space  5  is  reserved  for  the  number  N  of  |> '  s  used.  Since  the 
maximum  number  of  (j>  1  s  used  is  five,  0  <  N’<:  5,  Space  6  is  left  blank 
and  Spaces  7-12  are  reserved  for  the  weapon  in  millimeters. 

Spaces  13-18  are  reserved  for  the  projectile  weight  in  pounds. 

Spaces  19-24  are  reserved  for  the  chamber  volume  in  cubic 
inche  s . 
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Spaces  25-30  are  reserved  for  the  travel  in  inches. 

Spaces  31-36  are  reserved  for  the  maximum  pressure  in  psi. 

Spaces  37-42  are  reserved  for  the  charge  in  pounds. 

Spaces  43-48  are  reserved  for  the  muzzle  velocity  in  ft/sec. 

Spaces  49-54  are  reserved  for  the  web  in  inches. 

Spaces  55-62  are  reserved  for  the  burning  constant. 

Spaces  63-64,  65-66,  67-68,  69-70,  71-72  are  reserved  for  the 
value  s  of  <j>' s  . 

Card.  2  (This  card  is  only  required  if  no  propellant  code  - 

Spaces  2  and  3  above  -  is  given  or  is  given  as  "99.") 

Spaces  1-6  are  reserved  for  the  propellant  constant  a. 

Spaces  7-12  are  reserved  for  the  propellant  constant  a°. 

Spaces  13-18  are  reserved  for  the  propellant  density  in  lbs/in^. 

Spaces  19-24  are  reserved  for  the  propellant  co-volume  in  in^/lbs. 

Spaces  25-30  are  reserved  for  the  propellant  force  in  ft-lbs/lb. 

Most  of  these  values  must  be  punched  on  the  card  with’  a  decimal 
point.  The  exceptions  are: 

The  codes:  Spaces  1-4. 

The  burning  constant  which  may  be  given  to  eight  significant 
digits. 

The  <j>'s  less  than  unity.  The  values  for  <j>  are  0.05,  0.10, 

0.05,  ...,0.95,  1.00  and  are  punched  as  05,  10,  15,..., 

95,  1. 
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Output  Format 


The  output  will  list  all  input  data  as  well  as  the  calculated  loading 
densities,  tabular  values  and  velocities,  pressures,  web  or 

burning  constants  depending  on  the  problem  being  solved.  If  the 
problem  is  of  Type  1  or  2  and  the  (j>'s  given  as  input  do  not  result  in 
■optimum  efficiency,  consecutive  (j>'s  are  tried  in  an  effort  to  improve 
the  results. 

°  . 

<j>,  Pp  ,  Em,  ,  -2-  are  referred  to  as  "PHI,”  "PP(p,"  "XIM," 
"GAMMA"  and  "-2S"  respectively. 

An  example  of  this  program's  output  is  in  Appendix  B. 
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PART  II  -  NUMERICAL  CALCULATIONS 


Let  W  denote  the  weapon  (caliber)  in  millimeters.. 

Let  M  denote  the  projectile  weight  in  lbs. 

Let  V_  denote  the  chamber  volume  in  i'n^. 

Let  L  denote  the  length  of  travel  in  in  . 

Let  PMAX  denote  the  maximum  pressure  in  psi. 

Let  POPR  denote  the  operating  pressure  in  psi. 

Let  B  denote  the  burning  constant  in  (in/ sec) /psi. 

Let  C  denote  the  charge  in  lbs. 

Let  WEB  denote  the  web  in  inches. 

Let  Zb  denote  the  loading  density  in  gms/cc. 

Let  Vm  denote  the  velocity  in  ft/sec. 

The  following  are  constant  for  the  propellant.  Their  values  are 
given  in  Table  2  of  Appendix  A. 

Let  a  and  a°  denote  the  two  "propellant  constants"  in  in^/lb. 

Let  p  denote  the  propellant  density  in  lbs/in^. 

Let'll  denote  the  propellant  covolume  in  in^/lb. 

Let  F  denote  the  propellant  force  in  ft-lbs/lb. 


If  C  is  given  (Cases  2,  3  and  4),ZL(the  loading  density)  is  obtainable: 

Zb  =  C/Vc 

To  facilitate  the  calculations,  the  parameter  is  used  and  is  equal 
to  the  value:  a/  (  -  _1) 

From  |)  =  a/(^-  ~^)  >  the  equalities 

,11. 

■Hz*.  -  p)  -  a 

_1  _  a 

^  ~  p  ~  <t> 

]_  =]_  a 

^  p  <f) 


are  derived.  Therefore,  if  n  <J>' s  are  chosen,  ij)-^,  .  .  .  ,  (as  in  Case  1) 
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and  C  is  not  given,  each  of  the  n  s  are  calculable: 


+  7T ) 

p  <t>i 

and  hence  C^  =  i  *  Vc ,  i  =  1 ,  .  .  .  ,  n 
Let  A  denote  the  bore  area  of  the  projectile. 


A  =  ( 


A  )  -  ~rr 

2 


Let  X  =  (  Vc  +  -  A  •  L)  /  Vc 
Let  tm  =  (a^)/  (X  - 


Let  Ppi°  = 
PMAX  =1.15  POPR 


a°  PMAX 


M  +  3 


where 


Knowing  the  value  of  |>i  and  Ppi°,  the  values  of  Zi  and  '$■  i  for  each 
are  now  obtainable  from  Tables  4,  6,  8  and  10  in  Reference  1.  Tables 
4  and  6  give  -Zi  and  $~i,  respectively,  for  single-perforated  propellants 
Table  8  and  10  give  -Zi,  and.if'i,  respectively,  for  seven -perforated 
propellants. 


/  1.02(m  +  3.  lj 

Let  mi  =  32.  175 


The  velocity  Vmj  is  then  equal  to  the  value 


/i  *•  •  0.3  2CiF 

( 1  -  y  l  &m  )  — —  / 

0  •  3mj 


,  i  =  1 ,.  .  .  ,  n 


Consider  the  case  (Case  3)  when  C  and  Vm  are  given.  Here  Ci  =  C 
for. each  i  and  ij)i, P°,  , -Si ,  if"  i  and  .mi  all  have  the  single  values  cj>, 

Pp’  -Z ,  tf",  and  rr/  respectively. 


Therefore,  since  Vm 


=  l\j (1  -  fgm0'3) 


2  CF 
0 . 3m' 


„  2  0.3.  2CF 

Vm  -  (l-  )  0i3m, 


•,^m0.3  2£F/  = 


2CF 


0  ,  3m'  0 . 3m' 


Vm 


,  CF  ...  2.  0 .  15m 

'0 .  15m  '  ^  0.13 

£m  •  CF 


s 


(1  - 


0.  15m  Vm^  .  /  .g  q.  3 


CF 


)/  - 


Knowing  and  tj>,  Pp°  .  is  obtainable  from  Tables  6  or  10. 

From  Pp°  =  a°  PMAX 


M  +  i  ■* 


L  M  +  2  -i 


P? 


o. 


we  find  PMAX  =  ao 


£  1 

M  +  2 

c 

lM  +  3  J 


Knowing  Pp.  and  6,  Z  is  obtainable  from  Table  4  or  6 

t-  0. 3  c 

Let  U  m  =  S  • 

Let  (1  -  0.1485  Z)  /  &  _  =  Ts 

Let  (1  -  0. 242  Z)  /S'  =Tm 


When  a  single-perforated  propellant  is  used,  Ts  is  calculated -and 
the  point  of  optimum  efficiency  is  when  S  and  Ts  are  equal..  In  this  case 


'  WEB  =  (B/A) 


CFm  Z 


0.99 


However,  when  a  multiperforated  propellant  is  used,  Tm  is  calculated 
and  the  point  of  optimum  efficiency  is  when  S  and  Tm  are  equal.  In  this 
case  i - 


/  CF m 

WEB  =  (B/A)  Vi.  369  " 
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PART  III  -  PROGRAM  LOGIC 


Tables  8,  10,  4  and  6  of  Reference  1  appear  in  the  first  104  data 
cards.  As  I  =  1,  20  and  J  =  1,  47;  Cards  1-24  give  A(I,  J),  the  values 
in  Table  8.  Cards  25-48  give  C(I,  J),  the  values  in  Table  10.  Cards 
49-76  give  E(I,  J),  the  values  in  Table  4  and  Cards  77-104  give  F(I,  J) , 
the  values  in  Table  6.  I  =  Pp°/2,000  where  Pp°  =  2,000,  4,000,  6,000,. 
90,000,95,000,  100,000.  J  =  20\>  where  i  =  0 . 05  ,  0 .  10 , .  .  .  ,  0 . 95  ,  1 . 0 

Linear  interpolations  were  made  within  these  tables  when  2T  and 
-Z .  are  known  and  Pp°  is  to  be  found.  No  attempt  was  made  to  inter¬ 
polate  between  the  <j>'s;  instead,  the  tabular  closest  to  the  calculated  <j) 
has  been  selected. 
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PART  IV  --  EXAMPLES 

90mm,  Gun,  M41 


1.  Gun  constants 

Projectile  weight  M - 12.65  lbs. 

Chamber  volume  Vc  - - -  300  in.  ^ 

Propellant  Type - : - M17  (M.P.) 

Total  Travel  L - - - :  155  in. 

2.  Case  4 

Given:  Charge  C - . - 8,58  lbs. 

Maximum  Pressure  PMAX -  50,500  psi 

Velocity  Vm -  4,000  f/ s 

Web  - - - - -  0.052  in. 


Vc 


8.58 

300 


=  0.0286 


12.92 


12.92 


1 


1 

P 


1 


0.0286  "  0.0603 


34.965  -  16.5837 
12.92 


18. 38126 


0. 70289- ~  0.70 

(where  a  and  f3  are  found  in  Table  2  (9)  of  Appendix  A). 


Pp°  =  a°PMAX 


c 

f  8.58h 

M  +  *3  • 

=  (0.851)(50, 500) 

12.65+  -3  - 

c 

8.58 

_M  +  2  _ 

.12. 65  +  2 

=  (42975.5) 


15.51] 


.16. 94  J 


(42975. 5)(0. 915584)  =  39347.697 


(where  a°  is  also  found  in  Table  2). 
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=  1.763  by  Table  8  in  Reference  1. 

c  8.58 

/  _  1 . 02(M  +3.1  )  _  1.02  •  (12.  65  +_  3.1  ) 

m  32.174  '  32.174 

1.02(15.41777)  =  15.7261 

32.174  32.174 


=  0.48878 

90  2  2 

A=(^T4  *0.5)  (3.  1416)  =  (1.772)  ( 3.  14 16)  =  9.  8607 

'  /|  /  CFm'g, 

Since  Web  i  (b/a)  M/  1.369 

B  .  .  /!/  (8.  58)(364,000)  (0 . 48878)(  1 . 763) 

0.052=  9.8607  '  1.369 


, /2681157. 356 
■  0. 5127  =  B  *  '\l  1. 369 

r  1637. 424\ 
0.5118  =  B'V  1.17  ) 

0.598806  =  1637.424  B 
B  =  0.0003657  ' 


3.  Case  2 

Given:  Charge  C -  8.58  lbs; 

Maximum  Pressure  PMAX -  50,500  psi 


Zb  =  0.0286,’  <}>  =  0.70;  Pp°  =  39347.697 
^  =  1.410  by  Table  10  in  Hirschfelder 
A  =  9.  8607- 

X  =  Vc  +  A*L  =  300  +  (9. 8607)(155)  =  1828.4085 

Vc  300  300 

=  6.0947 
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_  (12. 92)(0.0286) 

°  m  6.0947  -  (29. 50)(0. 0286) 


0.3695  _  0.3695 

6.0947  -  0.8437  5.251 

S  =~t:  m°‘  3  =  0. 45  1 

Vm  =  y(l-  t€m0-3) 

11  1  0.  3m' 

=  j]/f\  -  (1.410)(0. 451)j  *  2(8. 58)  (364000)- 

J  (0. 3)(0. 48878) 

=  V(l  -  0.6359)*^— -4-6-— 

'  ''0.14669'' 


2274255.984 
0. 14669  ’ 


0.07037 


=  1508.4577  ■  =  3938.532 

0. 383 


1  -  0. 242  e  __  1  -  (0 , 242)(  1 , 763)  _  0,57393  =  Q  40? 

^  1.410  1.410  ' 


Since  S  =  0.451  >  0.407  =  Tm,  burning  is  taking  place  outside  the  weapon. 
To  attain  optimum  conditions,  new  and  smaller  <fc>'s  (and  hence  new  charges) 
•  may  be  chosen  in  an  effort  to  minimize  J  S  -  T-m|  .  This -was  done  on  the 
computer  and  the  results  are  in  Appendix  B  for  Cases  1  and  2. 

It  is  obvious  that  the  web,  calculated  with  B  =  0.000  3657,  will  be 
0,052  since  the  equation  for  web  was  used  in  Case  4  to  find  B.  ' 

4 .  Case  3 

Given:  Charge  C  - 

Velocity  Vm 


8.58  lbs. 
4,000  f/s 
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^  =  (1  l  0-  15m/  VmJ 2)  ^ 
C  F 


0.  3 


1  -  (0.  15)  (0.48878)(4000)2 
(8. 58)(:364,  000) 

0.451 

1173072 

1  -  3123120  =  1  -  0..  37561  =  0.62439  =  1.384 

0.451  0.451  0.451 


Since  <j>  =  0.70,  Pp°  =  41561.1367  by  Table  10  in  Reference  1. 
Therefore,  -Z  =  1.693 

c 


PMAX 


Pp° 


M  +  2 


£ 

M  +  3 


/  8.58 

/41561 . 1367  \  f 12.65  +  2 

'  (  0.851  /  1  8. 58 

' 12. 65  +  3 


/41561. 1367>i  r  16.  94  \ 
l  0.851  )  '  V  15.51  > 


=  (41561. 1367)(1. 2834) 


=  53340.792  psi 


Therefor  e  , 


Z  = 
Web 


1.693  by  Table  8  in  the  reference. 

0.0003657  .  J  (:8.58)(364000)(0.48878)(  1.693) 

9.8607  1.369 


0.0003657  '  a/258439597.29 

9.8607  1.17 


(0. 000 3657 )( 1607  6. 0566) 
11. 537 


0.05879 
11.  5  37 


=  0.005096 
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APPENDICES 


APPENDIX  A 


TABLES 


TABLE  1 


PROPELLANT  CODES 


Propellant 


Code 


Ml 

M2 

M5 

M6 

M9 

M 10 

M14 

M15 

M17 

T25 

T28 

T  34 

T  36 


1 

2 

5 

6 

9 

10 

14 

15 
17 
25 
28 
34 
36 


In  order  for  the  propellant  constants  to  be  read  in  as  part  of  the 
nput  (data  Card  2),  the  code  is  "99." 
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TABLE  2 


PROPELLANT  CONSTANTS 

In  the  following  table,  a  and  a°  represent  the  propellant  "a"  - 
constants;  P  represents  the  propellant  density  in  lbs/in^,^  represents 
propellant  covolume  in  in^/lbs.  ,  and  F  represents  the  propellant  force 
in  ft.  lb's  /lbs. 


Propellant 

a 

a° 

P 

x- 

.F 

Ml 

12.92 

1.015 

0.0567 

30.57 

305,000 

M2 

11.16 

.  0.743 

0.0597 

27.91 

360,000 

.  M5 

10.74 

0.725 

0.0596 

27.52 

355,000 

M6 

12.4/1 

0.938 

.  0.0571 

29.  92 

317,000 

M9  '  ' 

9.02 

0.566 

0.659 

25.97 

382,000 

M10 

11.15 

0.788 

0.0602' 

27.76 

339,000 

M14 

12.36  ■ 

0.906 

0.058 

29.54 

327,000 

M15 

14.  50 

•1.034 

0.06 

31.  17 

336,000 

M17 

12.92. 

0.851 

0.0603 

29.50 

364,000 

T25 

11.57 

0.786 

0.0585. 

28.66 

353,000 

T28 

11.68 

. 0.786 

0.0585 

28.77 

356,000 

.  T  34 

14.07 

1.007 

0.0596 

30 . 85 

335,000 

T  36 

12.59  ■. 

0.828 

0.06 

•29.26 

364,000  . 
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APPENDIX  B 


PROGRAM  OUTPUT  FOR  CASES  1-4 
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703  0.029  ■  8.58  .41561.1367  0.070  1.384  1.693  53340.792  0.451  0.426  0.050960869 
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APPENDIX  C 


FORTRAN  PROGRAM  FOR  INTERIOR  BALLISTICS 


FORTRAN  PROGRAM  F0R  HIRSCHFELDER  INTERIOR  BALLISTIC  CALCULATIONS 

DIMENSION  A(25,55),B(25),C(25,55),0(25) »  PH ( 1 5 ) , DELTA!  IS) ,  CH ( 1 5 ) , P0 
1P( 15) ,XIM( 15) , ZSC 15) .GAMMA (15) ,H(15),VM<15),S(15),T(15),WEB(15),E( 
225,5.5)  ,F  (25,55  )  ,  AB  (  50  )  ,  ABV  (  50  ) 

READ  INPUT  TAPE  2 , 2 , ( A ( 1 , J ) , J  =  l , 4 ) , ( A ( 2 , J ) , J  =  2 , 8 ) , ( A ( 3 , J  )  , J  =  3 , 1 3 ) , 
UA(4,  J )  ,  J  =  41l19  )j  (  A(5,  J  ),  J  =  5,24)  ,  (  A(6,  J),  J  =  5,3i  ),  (  A  (7,  J)  ,  J=6,38>  ,  (  A 
2 ( 8, J ) , J=7, 45 ) , (A (9, J ) , J=8,47) , ( A ( 10 , J ) , J=9 , 47 ) , ( A ( 1 1 , J ) , J = 10 , 47 ) , ( 

_ 3A(  12,  J  )  .  J=11.47)  .  (  A(  L3,.JJ,  J=12,47)  ,.(.A(  14,  J  ). ,  J  =  1 3 , 4  7  )  ,  (  A  (  1 5 ,  J  )  ,  J  =  1 4 

4,47>,(A(16,J),J=15,47),(A(17,J),J=16,47),(A(18,J),J=17,47),(A(19,J 
5) , J  =  18,47) , (At  20, J) , J=19,47) ,(B(I),I=1,5) 

2  FORMAT ( 24F 3.2/ (24F3.2 ) )  ' 

READ  INPUT  TAPE  2 , 3, ( C ( 1 , J ) , J  =  1 , 4 ) , ( C ( 2 , J ) , J  =  2 , 8 ) , ( C ( 3 , J  )  , J=3 , 1 3 ) , 

1 (C (4, J ) , J=4,19 ) , (C(5, J ) , J=5,24) , (C(6, J) , J=5, 31 ), (C(7, J ) , J  =  6, 38) , (C 

_ 2(8.J)«J  =  7.45).  (C(9,J).,J  =  8  ,47)  j  LC  ( .10.,.  J )  ,  J  =  9 ,47  ),(C1U  ,.J  )  ,3=10,47)  ,  ( 

3C< 12, J) , J=ll,47> , (C( 13, J) , J=12,47) , (C ( 14, J) , J=13, 47) ,  (C ( 15, J )  ,  J=14 

4.47)  , (C(16,J) , J=15,47) , (C( 17, J) , J=16,47) , (C(18,J) , J=17,47) , (C ( 19, J 
5 ) , J=18,47) , (C( 20, J) , J= 19,47) , (0< I ), 1=1,5) 

3  FORMAT ( 24F3- 2/( 24F3*  2)  ) 

READ  INPUT  TAPE  2 , 8 , ( E ( 1 , J ) , J= 1 , 7 ) , ( E ( 2 , J j , J=2 , 1 5 ) , ( E ( 3 , J ) , J= 3 , 25 ) 

_  _  1 ,  ( £ ( 4, J ) , J=4,3  5)  ,  <E(5,  J),U  =  4,46)  ,  (E(6,J),  J=5,47),  (E(7,-J)  ,  J  =  6,47)  ,  (._. 

2E ( 8 , J ) , J=7 ,47 ) ,(E(9,J) , J=8,47 ) , ( E ( 10 , J ) , J=9, 47 ) , ( E ( 1 1 , J ) , J=10 , 47 ) , 

3 ( £ ( 12 ,J) , J=ll, 471, (E { 13, J) , J= 12,47 ) , ( E t 14 , J ) , J=1 3 , 47 ) , ( E ( 15 , J ) , J=1 

44 . 47) , (  E  ( 16, J) , J= 15,47 ) ,(E(17,J),J=17,47),(E(I8,J) , J=18,47 ) , ( E( 19, 
5J  ).*J  =  19,47), (E  (20,J)  i  J  =  20,47) 

8  F0RMAT(24F3-2/ (24F3-2) ) 

_ READ  INPUT  TAPE  2,8, (F (1, J] , J  =  l, 7) , ( F I  2, J ) , J  =  2 , 15 ) , ( F ( 3 , J ) , J= 3, 25 ) 

1  ,  (  F  (  4 ,  J  )  ,  J.=  4, 3  5 )  ,  (FI  5,  J)  ,  J  =  4,46)  ,  (F  (6,  J)  ,  J  =  5,47)  ,  (  F  (  7  ,  J  )  ,  J=6 , 47  )  ,  ( 
2F(8,J) ,J  =  7,47)  , (F(9,J) ,J  =  8,47) , (F(10, J),J=9,47), ( F { 1 1 , J ) , J  =  10 ,47) , 

3 ( F (  12, J) , J  =  1 1 » 47) , (F( 13, J) ,J=12,47)  , (F( 14, J) , J  =  13,47) , (F( 15,J),J  =  1 

44.47)  , (F( 16, J) ,J=15,47 ) , (F ( 17, J) , J=17,47) , (F( 18, J ) ,J=18,47) , ( F( 19, 

5J) , J=19,47) , (F (20, J) , J=20,47)  , 

_ 9  FORMAT  (  24F3»2/  (24F3.2 )  )  . ......  . 

I  READ  INPUT  TAPE  2 , 4, K ASE , K0DE , M , N , W , P JW, VC , TR A V, P M AX , CH ( 1 ) , VM ( 1 ) , W 
1EB(1),Y,(PH(I)  ,1  =  1, N) 

4  F0RMAT( II, 12,2  II, l X , 3F 6 . 0 , F8 . 8 , 5F2 . 2 > 

P0PR=PMAX/1.15 

D0  500  1=1,20 

‘ _ A, Cl  »  46  )  =  A  (  1,47  ) 

7“  (  1,47)  =A  (1,50) 

Cl  I , 46 ) =C ( 1,47  ) 

C(  1,47) =C ( 1,50) 

Ell,  46  )  =£  (1,47)  3/ .13/64- 

E (  I  ,  47 ) =E (  1,50) 

_ _ F(  l  i_A6.1=F-l-I_.  4-7.1. .  I _  .  . . . . . . 

Ft  I ,47) =F ( I ,50 ) 

500  CONTINUE  Page  1 

IF(K0DE-99)6O»61,6O 

61  READ  INPUT  TAPE  2 , 62 , V , A0 , RH0 , ET A , U 

62  FORMAT ( 5F6.0) 

_ GO _ T.0...30Q _ _ _  ...  _  _ _ _ 

60  IF(K0DE-1)61,64,63 
64  V=12. 92 
A0= 1 . 01 5 
RH0=O . 0567 
ETA=30. 57 

..._ . . U=_305000. _  ..  _ _ _ _ _ _  _ _ _  _ _ - 
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63  IF  1K0DE-2X65,,6  6.  65 _ 

66  V  =  1 1  .  16 

A0  =  O  . 743  . 

RH0  =  0.0597 

ETA=.  2  7.91  . 

U=  360000.  . 

65-_  EE 1  K-0 D E -5  )_6  7._»_6_E* .6.7. . . . . ___  ■ 

68  V= 1 0 . 74 
A0  =  O. 725 

RH0  =  0.0596 
ETA  =  27.52 
U  =.  355000. 

67  IF (&0DE-6) 69, 70,  69 

70  V= 12 . 41 

A0  =  0.938 
RH0  =  0.0571 
ETA  =  29.92 
U  =  317000. 

69  1.F  ( K.0DE-9 )  7.1.,  7  2,  7.1 
72  V  =  9.02 

A0  .=  0.566 
RH0  =  0.6590- 
ETA  =  25.97 
U  =  382000. 

71  IF ( K0DE  -  10) 73,74,73 

74  V=  11.15 

A0  =  0.788 
RH0  =  0.0602 

ETA  =  27.76 

U  =  339000. 

-73  IF  (K0DE-141 7.5 , 76, 75 

76  V  =  12.36 
A0  =  0.906 
RH0  =  0.0582 
EtA  =  29.54 
U  =  327000. 

75  IF (K0DE-15) 77,  78,77 
78  V  =  14.50 

A0  =  1.034 
RH  0  =  0.0600 
ETA  =  31.17 
U  =  336000. 

77  IFIK0DE-17) 79, 80,79 

80  V=  12.92 
A0  =  0.851 
RH  0  =  0.0603 
ETA  =  29.50 
U  =  364000. 

— Z9  IF  lK0DE-^25.)  81,  82,81  ' _ _  _ _ 

82  V  =  11.57 

A0  =  0.786  Page  2 

RH0  =  0.0585 
ETA  =  28.66 
U  =  353000. 

81  IF  (K0DL-28J.83, 34,83 
84  V  =  11.68 
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_.  A0  =  0.780  _  _ _ 

RH0  =  0.0585 

ETA  =  28.77  . . . 

U  =  356000. 

83  I  F  (  K0DE-34  )  85.,  86  ,  85.  . .  . ..  _  ... 

86  V  =  14.07 

_ A0_...= _ 1  ..DQJ _ _ _ _ _ _ _ _ _ _ _ 

RH0  =  0.0596 
ETA  =  30.85 
U  =  335000. 

85  IF  (K0DE-36)3OO»88»3OO  . .  ....  _ 

88  V  =  12.59 

_ A0  =  Q-..82J3 _ : _ : _ _ 

RH0  =  0.060 
ETA  =.  29.26 
U  f  364500. 

300  STEV0  =  0.0 

IF1KASE-1 ) 100, 31, 100 

_ 1.Q0..LF.1KAS£-2JJLJQ2,  103^.1.0.2 _  _ _ _ ..... 

103  N=1 

DELTA ( 1 )=CH( 1 ) /VC 

PH!1>=V/! (1./ DELTA  ID )-( 1./RH0) ) 

G0  T 0  106 

600  1=2 

DELTA  I  I ) =CH LI ) /VC 

PH( I )=V/( ( l./DELTA! I ) ) - I 1./RH0) ) 

_ CH  (.2.)  =C.H  (  LI _  _ 

WEB ( 2 ) =WEB ( 1 ) 

_ V  M  (  21  =J/.MUL) _ _ _ _ _  _ _ _ _ _  .. 

G0  T 0  106 

105  DELTA ! 1 ) =  CH!  1)/VC 

PH( 1 )=V/ ! (l./DELTA! 1 ) )-( 1-/RH0) ) 

AREA  =  ( LW/Z5.4)**2)*(3. 1416/4.) 

X  =  (VC  +  ( AREA*TRAV) ) /VC 

_  XIM(L)  =  _(.V*DELTA(  1)  )  /  !X- LE1A*DELTA  (.1.1.11 _  _  _ _ 

S ( 1 )  =  XIM(1)**0.3* 

H(l)  =.  .0317*  (  PJW+!  .3226*CH(  1  )  )  ) 

GAMMA! 1)  =( 1 . — ( C . 15*H( 1 )*(VM( l)**2) )/ (CH! 1)*U) ) )/S(l) 

_ fl  =  LPH  ( J_1  *20  .  1+0.4 

L  =  Q 

_ I  7  =  0 _ _ _ _ _ _ _ 

D0  120  1=1,20 

_ IJFIL-.I  1120, 140,120  . . .  _.  _ 

140  12=1 

120  CONTINUE  Page  3 

IF ( 12 ) 122, 121 , 122 

_ 121.  WR  I  TE— 0.UTPUT  T  APJE  3 , 122., EASE ,RH.L1.) _ _ _ _ _ 

123  FORMAT! 1H1.10X  ,7HA  CASE  I3,63H  TYPE  PROBLEM  HAS  BEEN  REJECTED  HERE 
1  BECAUSE  CALCULATED  PHI  IS  F10.4.38H  FOR  WHICH  THE  TABLES  ARE  NOT 
2 ADEQUATE  ) 

_  GO  TO  40 

122  IF (M-l ) 108, 108 ,107 

_ .lQ7_qg  .1.0.9..  .  _L^.L,.4X _ _ _ _ 

ABV ( I )  =  ABSF ( C ( L  ,  I  )  - 
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109  CONTINUE  _ _ _ _ _ _ _ _ _ 

VAL  =  ABV! 1) 

D0  129  1=1,47 

VAL  =  M I N 1 F ( VA  L, ABV ( I ) ) 

129  CONTINUE  ..... 

D0  110  1=1,47 

IF  LAB  VTJLJ  -  VAU.  ,J_1Q  ,1.12  >1113.  - _ _ : _ _ _ .' 

112  K  =  I 
TK  =  K 

110  CONTINUE 

I  F  (  GAMMA  (  1 .)  .  -  C(L,K))113, 114,115 

113  R=(C(L,K) -GAMMA! T > ) / ( C ( L , K ) -C ( L , K+l ) ) + TK 

_ 1 _ G.0.  T0  116.  _ _ '  ■  . . . ■  _ _ _ _ _ 

114  R  =  TK 

30  T 0  116 

115  R= ( C ( L , K— 1 ) -GAMMA! 1 ) > / { C  t  L , K- 1 ) -C ( L , K ) l  +  TK-l . 

116  P0P(1)=R*2OOO. 

J  =  R 

ZStl)  =  A  (  L ,  J  )  — I  ABSF  (R-Z)*ABSf7a(L,  J1  —  A~(  L~Tj  +  lT7~) 

T(l)=.  (1.  — «242*ZS(1))/GAMMA(1) 

151  WEB!  1  )  =Y/ SQRTF  ( ( 1 . 369* AREA**2 ) / (CHI  1 ) *U*H ( 1 ) *ZS ( 1 )  ) ) 

PMAX=(P0P(1)/A0)*(1.+ICH(1)/(6.*PJW+2.*CH(1>  I  )  )  . 

G0  T 0  137 

. 108.00.130  1  =  1,47  _ _ _ _ _ _ _ _ _ _ _ _ 

ABV( I )=ABSF(F( L, I ) -GAMMA! 1 ) ) 

130  CONTINUE 
VAL=ABV ( 1 ) 

D0  131  1=1,47  . 

VAL  =  MIN1F(VAL*  ABV! I )  ) 

_ L3.1 _ C0-NT.I.N.UE _ _ _ _ _ 

D0  132  1=1,47 

IF(ABV(I)-VAL)  132,133,132  _ _ _  _ 

133  K= I 

TK  =  K  •  _  _  _ 

132  CONTINUE 

_ LTLL  GAMMA.! .  1  i  -iEjLL,  K1  11.50, .13  4 , 133 _ _ 

150  R=(F(L»K)-GAMMA(1) ) / ( F t L , K ) -F ( L , K  + 1 ) ) +  TK 

GO  T0  136  .  . . . . . .  . 

134  R=TK 

G0  T 0  136  .  .  . .  ....  .  _ _ _  _ _ _...  . 

135  R= (F (L.K-l ) -GAMMA! 1 ) ) / ( F ( L , K— 1 ) -F ( L , K ) )  +  TK-1 . 

..136  POP!  1  l=R»?00f>. _ _ 

J  =  R 

Z  =  J  .  .  ...  . . .  ...  ..  ..  .. 

ZSll)=E(L,J)-< ABSF(R-Z)*ABSF(E(L,J)-E(L,J+1} ) ) 

T(  1)  =  (  l.-.1485*ZSCin/GAMMA(l)  . 

161  WEB! 1 ) =Y/ SQRTF ( ( 0. 990* AREA**2 ) / ( CH( 1 ) *U*H( 1 )*ZS( 1 ) ) ) 

_ PMAXrr  (  POPiX]  /  A£  J.»  ( 1, ±.(.CH(1J  ZL6»*P  J H±2^*CH l  lJ  ll  ). .  . . . . 

137  WRITE  0UTPUT  TAPE  3 , 138  ,  KASE ,  K0DE  ,  M  ,  W  ,  PJW  ,  VC  ,  TRAV  ,  CH  (  1  )  ,.VM  ( 1 )  ,  ETA , 

1RH0,U,V,A0,Y,  AREA  ........ 

138  FORMAT ! 1H1 , 50X  ,9HI  N  P  U  T//5X,7HCASE  =  I 1/5X,  18HPR0PEL LANT  CODE  = 

1  I2/5X.21 HC0OE  FORM  FUNCTION  =  II, 6H  Z5X ,  16HWEAPCN  W  (MM)  =  F 

212.4/5X.22HPR0 J.WT.  PJW  (LBS.)  =  F 12 . 4/5X , 2 7HCHAMBFR  VOLUME  VC  (C. 
3 IN )  =  F12.4/5X.19HTRAVEL  TRAV.  (IN)  =  FJL2».4/5X ,^1RCHARGE  CH11.1  JIBS 
4)  =  F12.4/5X, 2 3H VELOCITY  VM(1)  (F/S)  =  F12.4.10H  /5X.32HP 
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5R0P.  C0V0LUME  £T  A  (C.  IN/LBS)  =  F  12.4./.5X,  3  1HPR0P.  DENSITY  RH0  (LBS/ 
6C.IN)  =  F12.4/5X.29HPR0P.  FORCE  F  (FT. LBS/LBS)  =  F 12 . 4/5X , 24HPR0P . 

. .  7.  C  0  N  S-T  AN  TS  A.A0,  =  F12.4, 5  X ,  FI  2, 4/5  X  ,18HB  URN. CONSTANT  B  =  F12.8/5X 

8* 12HB0RE  AREA  =  F12.4) 

_ WRITE” OUTPUT”  TAPE  3,8977pH( I ) .DELTA ( I ) ,CH( I ) , P0P (I ) ,XIM( I) .GAMMA  (  I 

_ 1  )  ■  ZS  t  I  )  .PMAX.S  (I  )  .  T  (  I  )  .  WEB  (  I  >  .  I  =  1  .  N  ) _ _ _ _ _ 

89  F0RMAT11X//5OX.L1H0  U  T  P  U  T//102H  PHI  L. DENSITY  CHARGE 

X  PPO  .  J<1M  GAMMA  ZS  PRESSURE  S  T 

2  WEB///I 1X.F8  .2, 1X,F8.3,1X,F10.2, IX.F12.4, 1X.F8 . 3.1X.F6.3, IX ,F6.  3 
___  3,1X,FL2,3.1X.F8.-3j_LX.EB.3j1X,FLZ.9)  ) 

G0  T 0  40 

,  153_.HR  I TE.  ..OUTPJJIXI.APX  3  »X5Ja_i  KASE_»-K2DE  iiM,W._PJR,.y£_,XRAV  ,  PM  AX,P0PR_,.£fcL(Xl  , 
1VM ( 1 ) « WEB ( L ) , ETA,RH0,U,V, A0, AREA 

L55  FORMAT ( IH1.S0X  «9HI  N  P  U  T//5X.7HCASE  =  I.1/5X  , 18HPR0 PELL ANT  CODE  = 
1  1 2/ 5  X  ,  2 1HC0DE  FORM  FUNCTION  =  II, 6H  /5X  ,  16HWEAP0N  W  (MM)  =  F 

212.4/5X, 22HP  R0  J . W  T .  PJW  (LBS.)  =  FI 2 . 4/5X , 27HCHAMBER  VOLUME  VC  (C. 

3  I N )  =  F12.4/5X .I9HTRAVEL  TRAV  (IN)  =  F 12 .4/5X , 26HMAX  PRESSURE  PMAX 
_ 4  (PSI)  =  F12.4/5X.32H0PERATING  PRESSURE  P OPR  (PSD  =  F 12 . 4/ 5X , 1 5HC 

5HARGE  (LBS)  =  FI 2 . 4/5X  ,  1 7H VEL 0C I T Y  (F/S)  =  F 1 2 .4/ 5X , 6HWEB  =  F12.4, 

654H  . .  /5X.32HP 

7R0P.  C0 VOLUME  ETA  (C. IN/LBS)  =  F 12. 4/5X, 3 1HPR0P.  DENSITY  RH0  (LBS/ 
8C.IN)  =  F12.4/5X.29HPR0P.  F0RCE  F  (FT. LBS/LBS)  =  F12-4/5X, 24HPR0P . 
9  CONSTANTS  A,A0,  =  F 12 .4, 5X , F 12. 4/5X , 1 2HB0RE.  AREA  =  F12.4) 

WRITE  OUTPUT  TAPE  3 , 700  ,PH(2),DELTA  <  2  ),CH(1  )  ,  P0P"(  2  )~,  X  I M  (2  )  .GAMMA  (2 
1).ZS(2),VM(1),S(2).T(2).Y 

700  F0RMAU 1X//5OX.11H0  U  T  P  U  T//102H  PHI  L. DENSITY  CHARGE 

1  PP0  X I M  GAMMA  ZS  VELOCITY  S  T 

2  B///(IX,F8.2,1X,F8.3,IX,F10.2,1X,F12.4,1X,F8.3,1X,F6.3,1X,F6.3 

_ 3. 1X,F12.3,LX,F8.3,1X,F8.3, 1X.F12.9) )  . .  . 

WRITE  OUTPUT  TAPE  3,156 

156  FORMAT ( 1X//30X .57HTHE  POINT  OF  OPTIMUM  EFFICIENCY  IS  WHEN  S  AND  T 
I  ARE  EQUAL) 

104  GO  TO  40 
31  DO  5  I  =  l.N 

DELTA  (  I  ).  =.  1./ (t  V/PH(I)  )  +  (  1./RH0)  ) 

CH ( 1 ) =VC*OEL  T  A ( I ) 

106  P0P(I)=  PMAX*A0* { l PJW  +  CH(  I  )/3.  )/ (PJW  +CH(I)/2.)) 

AREA  =(  (W/25. 4  )**2)»( 3.1416/4. ) 

X  =(VCMAREA*TRAV)  1/VC 

XIM( I )  =  (V»DELTA(  I  )  )/(X-(ETA*DELTA( I ) ) ) 

. . JI=JLP_H1X1*JLDQ  .1/5* _ _ . _ . . .  ..  ..  . . 

Q  =  Q+0 .4 

R=  P0PI I J/2000 . 

R  =  R+(  l./( 10. **7 ) ) 

L  =  Q 
J  =  R 

_ Z=J _ _ _ _ _ _  _  _  Page  5  _ 

I F ( M- 1)11,11,10 

11  ZS ( I ) =E ( L »  J ) - ( ABSF ( 1  POP ( I  1 /2 000. ) -Z ) * ABSF ( E 1 L , J ) -E l L ,  J+l) ) ) 

GAMMA  (  I  )  =F ( L , J  )- ( ABSF (  ( POP ( I ) /200  0 . )- Z )* ABSF ( F ( L , J )  — F ( L , J+l )  )  ) 

H ( I  )  =  . 031 7* ( P J  W  + l .3226*CH( I ) ) ) 

VM(  I )=SQRTF (  ( 1  .-GAMMA (  I )*(  X I  M  (  I  ) **0 . 3 ))*(( 2 . *CH ( I  )  *U )/( 0 . 3*H (  I  )))  ) 

_ S.L1 )  =.X  1  Ml IJ  »  *0  . 3  _ _ _  _ _ . 

T( I  )=( l.-< .148  5*ZS( I  )  )  ) /GAMMA ( I  ) 
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WFRI  I  l=Y/SORTF((Q.  990*  ARF  A*  *2  1  /  (  CH  (  I  )*U«H(  H«ZS(  I  )  )  )  . . 

I F { KASE-4 )  5,162,162 

162  Y  =  WEB(1)*SQRTF(( 0* 99* ARE A* *2 ) / ( CH ( 1 ) *U*H ( I  ) *ZS (  I  )  )  ) 

G0  T 0  153 

10  ZS(  I  )=A(L,J)-(  ABSF(  (P0P( I ) /2000 . ) -Z ) * ABSF ! A ( L  ,J)-A(L,J+i))> 

GAMMA  (  I)=C(L,  J  )-{  ABSF  (  (  P0P  (  I  )  /2000.  ) -L  )  *  ABSF  (  C  (  L  ,  J  )— C  (  L  ,  J  *■  1 )  )  ) 

_  _H(  Ii=.A31  7*  ( PJ.W  .+..  1-.3226*CH  U  J  J.)  . . . _ . . _ . . . 

VM( I ) =SQR  TF ( ( 1 .-GAMMA ( I ) *( XIM( I ) **0.3 )  )*(  ( 2 . *CH ( I  ) *U I / ( 0 . 3*H (  I) ) )  ) 
S(I)=XIM(1)**0.3 

T  t I )  =  ( 1 . - ( .  242  *Z  S (  I  )  )  ) /GAMMA (  I ) 

WEB( I )=Y/SQRTF ( ( l .369*AREA**2)/( CH( l ) *U*H(  l  )*ZS( I  )  )  ) 

I F ( KASE-4 )  5,152,152 

152  Y  ....=  WEB  I  J.  1  *SQR  TF  LLL  36  9*  ARE  A*  *2 )  /  IC.HI  1 )  *U*H(  I  )  *ZS  ( I )  1  ) 

G0  T 0  153 
5  CONTINUE 

IF(STEV0)2O1,2  00 ,201 

200  WRITE  OUTPUT  TAPE  3 , 6 , KASE , K0DE , M , W , P J W , VC , TR A V , P MAX , P0PR , ET A , RH0 , 
1U, V,A0, Y.AREA 

-6.  .FORMAT  L1H1,5QX_,9HI„N  „E_  U  JO./5X,.7HCASE  =_. LL/3.X.,  18HPR0  PELLANT  CODE  = 

1  I2/5X.21HC0DE  F0RM  FUNCTION  =  II, 6H  /5X , 16HWEAP0N  W  (MM)  =  F 

212.4/5X.22HPR0J.WT.  PJW  (LBS.)  =  F 1 2 . 4/5X , 27HCHAM BER  VOLUME  VC  (C. 
3TN )  =  F12.4/5X  .19HTRAVEL  TRA V  (IN)  =  F 12 . 4/ 5X , 26HM AX  PRESSURE  PMAX 
4  (PSI)  =  F12.4/5X.32H0PERATING  PRESSURE  P0PR  (PSI)  =  F 12 . 4/5X , 32HP 
5R0P .  C0V0LUME  ETA  (C. IN/LBS)  =  F 1 2 . 4/ 5X , 3 1 HPR0P .  DENSITY  RH0  (LBS/ 
6C.IN)  =  F12.4/5X.29HPR0P.  FORCE  F  1FT. LBS/LBS)  =  F 12 . 4/5X , 24HPR0P . 

7  CONSTANTS  A, AO,  =  F12. 4, 5X,F12.4/5X,18HBURN. CONSTANT  B  =  F12.8/5X 
8 , 12HB0RE  AREA  =  F12.4) 

WRITE  OUTPUT  TAPE  3, 20  ,  (  PH (  I  )  ,  DELTA ( I ) ,CH ( I ) , P0P (I) , X  I M ( I ) , GAMMA ( I 
1),Z_S(I),VM(I),S(I),T(I)  ,  WEB  (  I  ),  I  =  1,N) 

20  FORMAT ( IX/ /50X ,1 1H0  U  T  P  U  T//102H  PHI  L. DENSITY  CHARGE 

.1  PP2  _  XIM  GAMMA  ZS  VELOCITY  S  T . 

2  WEB/// ( 1X»F8.2»  1X»F8.3»1X,F10.2»1X»  F12.4, IX»F8.3,1X,F6.3,1X,F6.3 
3,1X,F12.3,1X,F8.3,1X,F8.3,1X,F12.9)) 

WRITE  0UTPUT  TAPE  3,7 

7  FORMAT! 1X//30X .57HTHE  POINT  0F  OP  1 1  MUM  EFFICIENCY  IS  WHEN  S  AND  T 
1ARE  EQUAL) 

IF.  I STE  VO  )  40, 3.2  ,4  0  . . . '  ...  ..  ..  ...  . . 

32  N4=2 
N5  =  1 

D0  30  I  =  1 ,  N 
I F ( S ( I ) - T ( I ) ) 29, 27,25 
25  AB( I ) =100. 

N5=2  _ _ _  _ _ _ _ _ _ _ 

G0  T 0  30 

27  WRITE  OUTPUT  TAPE  3 , 2 8 , PH ( I ) , S (I) 

28  FORMAT ( 1X//30X ,6HPHI  =  F12.3.15H  SINCE  S  =  T  =  F12.4) 

N4  =  1 

GO  TO  30 

29  AB(I)  =  ABSFJSII)  -  T.LI)  L  Page  6 

VALUE  =  AB ( 1  ) 

30  CONTINUE 
G0  T  0 ( 40 , 35 )  , N  4 

35  DO  23  1  =  1, N 

VALUE  =  MIN1F( VALUE  ,AB(  I  )  ) 

23  CONTINUE 
K=  10 
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D0  24  J=ljN _  .  .,  _  . 

IF (AB(  I  > -VALUE  >24,26,24 
26  K=  I 
24  CONTINUE 

I F  (  K-  IQ  )  950 ,95  1,  9.5 1 

951  WRITE  0UTPUT  TAPE  3 , 95 2 , PH (1)  ,  S (  1  )  , T (  1) 

952  F  0RMAT  (  1 X/  /  2QX  .22HTRU.E  PHI  IS  LESS.  THAN  F12.4.11H  WHERE  S  =  F12.4, 
1,9H  AND  T  =  F12.4//) 

961  PHI  1  )  =  P.H(l)  -  0.05 

IF (PH ( 1 )-0. 05) 40, 970,970 
.970  M=  1 

STEV0=2.O 

G0  T 0  31 _  _ _ _ 

950  IF ( STEV0) 900 , 9  05  » 900 
900  IF(K-2)201, 201.902 

905  WRITE  OUTPUT  TAPE  3 , 33 , PH ( K ) , S ( K  )  , T ( K ) 

33  FORMAT! 1X//20X .25HTRUE  PHI  IS  GREATER  THAN  F12.4.11H  WHERE  S  =  F12 
1.4.9H  ANO  T  =  F12.4//I 

_.9.02 I-FJLRH  (  K  )  -1  .  Q14_1..4.Q^4H . . . . . . .  .  _ _ 

41  IF(PH(K)-PH< 1)  >40,42.43 

42  IFIPHI l)-0.95) 50,51.51 

51  PH ( 1  )  =  1.0 
N=  1 

G0  T  0  48 

_.  5(1. J.  FCPHU1-Q . 901. 53, 52,52.  . .  . ... 

52  PHI  1 1=0.95 
PH (21=1.0 
N  =  2 

G0  T 0  48 

53  PH  I  1 )  =  PH ( 1 )  +  0.05 

PJil.2  )  =PH  ( 1 )  _t  C.  Q5 
PH ( 3 )  =  PHI2)  +0.05 
N=  3 

G0  T 0  48 

43  IF.(PH(K)-0. 95)  44,45,45 

45  PHI  1)  =  PHIK)  +  0.05 
N=  1 

G0  T 0  48 

44  IFIPHIK)  -  0.90)46,47,47 

47  PHI  1)  =  PHIK)  +  0.05 

PH  I  2 )  =  PH  1 1 )  +  0.1 

N  =  2 

_ G2  T0  48 . . . . . .  . .  ...  ........ 

46  PHI  1)  =  PHIK)  +  0.05 

PHI  2 )  =  PH ( 1 )  +  0.05 

PH  I  3 )  =  PH  I  2 )  +  0.05 

N=  3  Page  7 

48  STEV0  =  1.0 

_ S.0...I0_3.1  _ _ _ _ 

201  WRITE  OUTPUT  TAPE  3, 2 02 , I  PH  I  I  ) , DELT A  I  I  )  , CH I  I  )  ,  P0P I  I  )  ,  X  I  M I  I  )  , G AMM A  I 
II) ,2Sl I ) ,VM( I )  ,S(  I  )  ,T( I ) .WEB  I  I ) , I=1,N) 

202  FORMAT! 1X//(1X,F3.2,1X,F8.3,1X,F10.2,1X,F12.4,1X,F8.3,IX,F6.3,1X,F 
16.3, 1X.F12.3, 1X.F8.3, 1X.F8.3, 1X.F12.9) ) 

IF  I STEV0-1 . ) 40 ,92  0,960 
960  IF  I  SI  LJ.-TI  1.)  )  40,  27, 961 
920  D0  921  1=1, N 
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I F ( S 1 1 J -T  m I  921, 921 , 922 
922  N5=2 
921  CONTINUE 

IF ( N5-2) 32 » 40 , 40 
40  G0  T 0  1 

END (1,1, 0,0, 0,0, 1,1.0, 0,0, 0,0. 0,0) 
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